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Abstract. The origin of the azimuthal anisotropy in particle yields at high pr (pt > 5 GeV/c) in RHIC 
collisions remains an intriguing puzzle. Traditional flow and parton energy loss models have failed to 
completely explain the large V2 observed at high pr - Measurement of this parameter at high pr will help 
to gain an understanding of the interplay between flow, recombination and energy loss, and the role they 
play in the transition from soft to hard physics. Neutral mesons measured in the PHENIX experiment 
provide an ideal observable for such studies. We present recent measurements of 7r° yields with respect to 
the reaction plane, and discuss the impact current models have on our understanding of these mechanisms. 
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1 Introduction 

Two of the greatest mysteries that have arisen from the 
RHIC physics program are the source of the apparent flat- 
ness of the high pr (> 5 GeV/c) suppression of Raa E3 
and the source of non-zero V2 at high px |Ej- The exis- 
tence of intermediate to high pr «2 was suggested early in 
the RHIC program 0, and has been the subject of many 
theoretical treatments (see 0]G3 for some additional ex- 
amples). Traditional flow and parton energy loss pictures 
have failed to describe the magnitude of this anisotropy. 
Measurement of the azimuthal asymmetry V2 at high px 
will shed light on the contributions from flow, recombina- 
tion, and energy loss, as well as the transition from soft 
to hard production mechanisms. 



2 Measuring v 2 and n° yields in PHENIX 

The orientation of the reaction plane is measured event- 
by-event using the set of two PHENIX Beam-Beam Coun- 
ters (BBCs), which reside at the region 3 < \r]\ < 4. Each 
detector is an array of 64 hexagonal, close-packed quartz 
Cherenkov counters, located 150 cm from the interaction 
point. The charge measured by each counter is propor- 
tional (on average) to the multiplicity of particles hitting 
it. The reaction plane angle ^rp is determined from the 
value of (cos 20). Because the two BBCs provide indepen- 
dent measurements oi^Ppp, we can estimate the resolution 
of the combined measurement via standard techniques 6 . 

For measuring photons and 7r°s, we use the Electo- 
magnetic Calorimeter (EMCal) 7 . Candidate clusters are 
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required to pass 7 identification cuts, and rriinv distribu- 
tions are formed from pairs of these clusters. The result- 
ing yields are binned in angle with respect to the reaction 
plane (A<f> — <p — &rp). A similarly binned mixed event 
background is then subtracted. The counts in the remain- 
ing peak centered on the ir mass are integrated in a ±2er 
window (where a is the width of a Gaussian fit to the 
peak). Six bins in A<j> are used in the interval [0 — vr/2]. 

To measure i>2, we fit the raw (uncorrected) A<f> distri- 
bution Y(A(f>) as 



Y raw (A(f>) oc l + 2v 



ran; 
2 



cos(2Z\0). 



(1) 



The resulting V2 parameter needs to be corrected for the 
reaction plane measurement resolution, hence the desig- 
nation w™" 1 . The resolution app is determined for each 
centrality bin, and leads to the corrected value V2° rr = 
v rau ' J aRp . The yields as a function of A<p can then be 
corrected with a factor 



Y(A4>) =Y raw {Acp) x 



1 



2v corr 



cos 2A(p 



2v r 2 aw cos 2A(f> 



(2) 



3 Results and Discussion 

To obtain Raa{A4>), we exploit the fact that the ratio of 
the yield at a given A<f> to the inclusive yield is equivalent 
to the ratio of the angle-dependent Raa to the inclusive 
Raa- Thus multiplying these relative yields by an inclu- 
sive measured Raa, we have: 



R AA {A<f>)=Y{A<f>)/Y x Raa 



(3) 



The Raa(A<P,pt) as a function of (N part ) is shown 
in Figure H We note that there appears to be a slightly 
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Fig. 1. Raa(pt) as a function of {N pa rt) for Au+Au 200 
GeV/c collisions; Each curve corresponds to a different Acj> bin 
and the panels are for different pr bins. The grey bands indi- 
cate the systematic error due to the reaction plane resolution 
correction. 
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Fig. 2. 7T° V2(pr) for Au+Au 200 GeV/c collisions, with each 
panel corresponding to a centrality bin. The grey bands indi- 
cate the systematic error due to the reaction plane resolution 
correction. 



different centrality dependence, according which A<f> bin 
in which the Raa is being measured. This feature is em- 
phasized by plotting the data on a semi-log scale, showing 
that the Raa behaves differently in different Acj> bins. 

The resulting ir° v-i is shown in Figure [2] For the first 
time we observe V2 up to 10 GeV/c. While the value of 
the V2 decreases beyond intermediate pr, it nonetheless 
shows a substantial and perhaps constant value out to the 
highest measured transverse momenta. 

To gain insight into the i>2 mechanisms at work at 
high pt, we turn to models. We compare the ir° V2 to 



two models, a calculation done by Turbide et al. [Hj (us- 
ing an Arnold-Moore-Yaffe (AMY) formalism 9 ) and the 
Molnar Parton Cascade (MPC) model ^01- Figure|3]shows 
calculations from these models, plotted alongside data for 
similar centralities. The AMY calculation contains energy 
loss mechanisms only, and we see that the data appear 
to decrease to a value at high px that is consistent with 
this model; the level of agreement is most striking in the 
20-40% bin. 

The MPC model has a number of mechanisms, includ- 
ing corona effects, energy loss, and the ability to boost 
lower partons to higher px (a unique feature). The cal- 
culation shown in FigureOHdoes a better job of reproducing 
the overall shape of the , though it is systematically low. 
It is important to note that this calculation is done for one 
set of parameters, so it should be very interesting to see if 
the MPC can better reproduce the data for a different set 
of parameters (the opacity of the medium, for example) . 

The prevailing thought is that the high pt behavior of 
the x>2 is due to energy loss mechanisms. If this is true, 
the Raa should be sensitive to the geometry of the col- 
lision. To test this behavior, we seek to combine the two 
traditional geometric parameters (centrality, or collision 
overlap, and angle of emission) into a single parameter, a 
quantity which we will refer to as "p L dV . Details of the 
calculation are described in as well as below. 

The Guylassy-Levai-Vitev (GLV) formalism can be used 
to calculate jet energy loss for a set of scattering cen- 
ters {xi}, where Xi = (ti, rj). In practice, an average over 
these centers is performed. As shown in Figure 01 if a 
static uniform color charge density within some region 
(p(x) — po and zero outside the region) is assumed, the 
resulting energy loss is AEqcd <x po-^max- More realis- 
tically, if the density seen by the particle changes along 
the path, we have AEqcd oc J q p{L)LdL (which reduces 
to the quadradic L dependence for a constant density). 
Application of this to a ID Bjorken expansion, with 



P(r, T ) = p(r,T ) — 

T 



(4) 



and given a jet trajectory r(r) = ro + v(r — To) (assuming 
v ~ c = 1 for the jet), we have 



L{t) = |r(r) i"o | — t To 



Therefore we have 



AE. 



QCD OC 



/' 

Jo 



p(r, r)LdL 

p(r,T )T L 
tq+L 



dL 



(5) 

(6) 
(7) 



This effective energy loss is calculated from the parton- 
density weighted average of the length from hard-scatter- 
ing origin to edge of an ellipse. Additionally, we perform 
a Glauber Monte Carlo sampling of starting points to ac- 
count for fluctuations in the location of the hard-scattering 
origin of the particles' paths within the region of overlap 
between the colliding nuclei. The crucial feature of p L dL 
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Fig. 3. Comparison of 7r° V2 with models. The top three panels show the AMY calculation with data for three centralities 
The bottom two panels compare two centralities with the b — 8 fm calculation of the MPC model |10) . 




Fig. 4. Schematic of the calculation of J p L dL. See text for 
details. 



is that it is proportional to the energy loss sustained by 
the parton as it traverses the medium. 

The resulting dependence of Raa for all centralities 
and angles on p L dL is shown in Figure [3] If the observed 
Raa arose from only geometric effects, we would expect 
the data to exhibit a universal dependence on p L dL. For 
low px, this is clearly not the case; something more than 
just energy loss is taking place there. However, when the 
Pt reaches 7 GeV/c and above, the Raa data do indeed 
appear to have a dependence on a single p L dL curve. 
This apparent scaling strongly suggests that the dominant 
effect on Raa at high-p-r is energy loss. These data can 
help to constrain energy loss models, and perhaps help to 
understand the nature of that energy loss (is it radiative, 
collisional, or some combination of both?). 



4 Conclusions 

We have presented the first measurement of high pt vi for 
7r°s. It is now clear that the v-i at high pt does decrease 
but to a non-zero value. Comparison of vi with models 
suggest that the dominant mechanism at work at high 
Pt is energy loss. In addition, we have presented the first 
measurement of ir° Raa as a function of angle with respect 
to the reaction plane. When the Raa data are examined as 
a function of an effective path length through the medium, 
the scaling that arises at high pr also argues for energy 
loss as the dominant mechanism at work. 
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Fig. 5. _Raa(^0,Pt) vs. p L dL. The panels correspond to different pr ranges. The solid circles are the most peripheral events, 
while the solid stars are the most central events. 



